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ABSTRACT 

We have measured the clustering of 0.30 < z < 0.90 red galaxies and con- 
strained models of the evolution of large-scale structure using the initial 1.2D° 
data release of the NOAO Deep Wide-Field Survey (NDWFS). The area and 
B W RI passbands of the NDWFS allow samples of > 10 3 galaxies to be selected 
as a function of spectral type, absolute magnitude, and photometric redshift. 
Spectral synthesis models can be used to predict the colors and luminosities of 
a galaxy population as a function of redshift. We have used PEGASE2 mod- 
els, with exponentially declining star formation rates, to estimate the observed 
colors and luminosity evolution of galaxies and to connect, as an evolutionary 
sequence, related populations of galaxies at different redshifts. A red galaxy sam- 
ple, with present-day rest-frame Vega colors of Bw — R > 1.44, was chosen to 
allow comparisons with the 2dF Galaxy Redshift Survey and Sloan Digital Sky 
Survey. We find the spatial clustering of red galaxies to be a strong function 
of luminosity, with r increasing from 4.4 ± 0.4/i _1 Mpc at Mr — 5\ogh ~ —20.0 
to 11.2 ± 1.0/i _1 Mpc at Mr — b\ogh fa —22.0. Clustering evolution measure- 
ments using samples where the rest-frame selection criteria vary with redshift, 
including all deep single-band magnitude limited samples, are biased due to the 
correlation of clustering with rest-frame color and luminosity. The clustering of 
—21.5 < Mr — 5\ogh < —20.5, Bw — R > 1.44 galaxies exhibits no significant 
evolution over the redshift range observed with r = 6.3 ± 0.5/i _1 Mpc in comov- 
ing coordinates. This is consistent with recent ACDM models where the bias of 
L* galaxies undergoes rapid evolution and r evolves very slowly at z < 2. 

Subject headings: cosmology: observations — large-scale structure of universe - 
galaxies: elliptical and lenticular, cD 
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1. Introduction 

The strength and evolution of the spatial clustering of galaxies is a function of cosmology, 
galaxy mass, and galaxy formation scenarios in Cold Dark Matter models (e.g., Peacock 
1997; Kauffmann et al. 1999; Benson et al. 2001; Somerville et al. 2001). The predicted 
linear and quasi-linear growth of density perturbations results in the spatial clustering of 
dark matter undergoing rapid evolution at z < 2 (e.g., Peacock 1997). However, CDM 
semi-analytic models and simulations predict the distribution of M R < —20 galaxies to be 
highly biased with respect to the underlying dark matter distribution (e.g., Cole & Kaiser 
1989; Kauffmann et al. 1999; Benson et al. 2001; Somerville et al. 2001). CDM models for 
an Q m = 0.3, A = 0.7 cosmology (ACDM) predict little evolution of the early-type galaxy 
spatial correlation function (measured in comoving coordinates) at z < 2. 

At low redshift, strong constraints on galaxy clustering as a function of spectral type and 
absolute magnitude are provided by the 2dF Galaxy Redshift Survey (2dFGRS) and Sloan 
Digital Sky Survey (SDSS) wide-field spectroscopic surveys (Norberg et al. 2002; Zehavi et al. 
2002). At z > 0.3, the comoving volume and galaxy counts of spectroscopic surveys are orders 
of magnitude less than the 2dFGRS and SDSS and the resulting estimates of the spatial 
correlation function suffer from small number statistics and biases from unrepresentative 
sample volumes (Hogg, Cohen & Blandford 2000; Giavalisco & Dickinson 2001). Only wide- 
field surveys, with large samples (> 10 3 ) of galaxies and volumes of > 10 6 /i~ 3 Mpc 3 , are able 
to provide robust estimates of r (the spatial scale where the correlation function equals 1) 
with errors of < 10%. 

It is possible to obtain constraints on the clustering of high-redshift galaxies using the 
two-point angular correlation function and deep wide-field imaging. While survey volumes 
at high redshift have increased with the advent of large-format CCDs, significant issues 
remain in the interpretation of angular correlation functions derived from these samples. 
Estimates of the spatial correlation function derived from imaging rely on models of the 
redshift distribution which differ significantly from each other. In addition, the distribution 
of galaxy types and luminosities can be a strong function of apparent magnitude and redshift. 
As the clustering of galaxies is a function of spectral type and luminosity (e.g., Davis & Geller 
1976; Loveday et al. 1995; Norberg et al. 2002; Zehavi et al. 2002), the apparent evolution 
of the correlation function can be dominated by selection effects rather than evolution of 
large-scale structure (Efstathiou et al. 1991). 

The advent of reliable photometric redshifts allows improved constraints on the redshift 
distribution of faint galaxies in imaging surveys (e.g., Brunner, Szalay & Connolly 2000; 
Brown, Boyle & Webster 2001; Teplitz et al. 2001; Firth et al. 2002). However, as discussed 
in §5, the resulting model redshift distributions and spatial correlation functions strongly 
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depend on the uncertainties of photometric redshifts. If rest-frame color criteria are also 
applied, samples of galaxies containing a comparable range of spectral types can be selected 
over a broad range of redshifts (Brown et al. 2001, Firth et al. 2002; Wilson 2003), thus 
reducing the selection effect which dominates clustering evolution estimates derived from 
single-band imaging. Absolute magnitude selection criteria can also be used, though this 
requires accurate photometric redshifts. 

To measure the evolution of galaxy clustering, we have used optical images from the ND- 
WFS (Jannuzi & Dey 1999, Jannuzi et al., in preparation) to select comparable samples of red 
galaxies at multiple epochs. We have used photometric redshifts and the PEGASE2 galaxy 
spectral evolution models (Fioc & Rocca-Volmerange 1997) with exponentially declining star 
formation rates to estimate the redshifts, luminosities, and rest-frame colors of all galaxies 
detected in the B W) R, and /-bands. In addition, we used the best-fit PEGASE2 models 
to estimate the spectral evolution and the present-day rest-frame colors and luminosities. 
We selected a red galaxy sample, with present-day rest-frame colors of Bw — R > 1.44, to 
allow direct comparison of the NDWFS clustering measurements with the low redshift early- 
type galaxy clustering measurements from the 2dFGRS. We also selected a Bw — R > 1-77 
subsample to allow a measurement of galaxy clustering as a function rest-frame color. 

The outline of the paper is as follows. We provide a brief description of the NDWFS 
imaging data and catalogs in §2. We describe the estimation of photometric redshifts and 
the galaxy spectral evolution models in §3. In §4, we discuss the motivation for studying 
red galaxies and the selection of these objects. We discuss the measurement of the angular 
and spatial two-point correlation functions in §5. In §6, we present the clustering of 0.30 < 
z < 0.90 red galaxies as a function of absolute magnitude and redshift. We discuss the 
implications of our results in §7 and summarize the paper in §8. 

2. The NO AO Deep Wide-Field Survey 

The NDWFS is a multiband (B w , R, I , J, H, K) survey of two « 9D° high Galactic 
latitude fields with the CTIO 4m, KPNO 4m, and KPNO 2.1m telescopes (Jannuzi & Dey 
1999). A thorough description of the observing strategy and data reduction will be provided 
by Jannuzi et al. (in preparation). This paper utilizes the first 1.2D° B W RI public data 
release of four adjacent KPNO 4m MOSAIC subfields in the Bootes field. The coordinates, 
depth, and image quality of the four subfields are provided in Table 1. This dataset, along 
with explanatory material, is available from the NOAO Science Archive on the World Wide 
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Web 1 . 

We generated object catalogs using SExtractor 2.2.2 (Bertin & Arnouts 1996) run in 
single-image mode with the minimum detection area, convolution filter, and signal above 
sky threshold optimized to provide the deepest catalogs possible as a function of the seeing. 
Detections in the different bands were then matched using criteria based on the distance 
between the image centroids. At faint magnitudes (R > 22), images with centroids within 1" 
of each other in different bands were matched. At bright magnitudes, images were matched 
if the probability of the centroids being within a given distance of each other by random 
chance was < 0.005. This probability was determined using the number of objects as a 
function of magnitude per unit area measured from the NDWFS. Where multiple matches 
were found, the closest centroids were matched. In practice, such criteria work well for both 
faint objects and bright objects including saturated stars and well resolved galaxies. 

We have used SExtractor MAG_AUTO magnitudes (Bertin & Arnouts 1996), which 
are similar to Kron total magnitudes (Kron 1980), due to their low uncertainties and small 
systematic errors at faint magnitudes. The la uncertainty of the MAG_AUTO photometry 
decreases from ~ 0.3 magnitudes at the 50% completeness limit to ~ 0.1 magnitudes at 2 
magnitudes brighter than the completeness limit. For unsaturated objects, the MAG_AUTO 
photometry has systematic errors which are significantly less than the la uncertainties. The 
measured correlation functions do not strongly depend on the technique used to measure 
object fluxes. Correlation functions using galaxies selected with MAG_AUTO photometry 
and 5" aperture photometry differed by < la. Throughout the remainder of the paper, all the 
results and conclusions are derived from NDWFS samples using MAG_AUTO photometry. 

We determined the completeness as a function of magnitude by adding artificial stellar 
objects to copies of the data and recovering them with SExtractor. The 50% completeness 
limits vary between the 4 subfields in the ranges of 26.2 < B\y < 26.8, 24.8 < R < 25.3, 
and 23.6 < I < 25.6. We used SExtractor's star-galaxy classifier to remove objects from 
the galaxy catalog which had a stellarity of > 0.7 in 2 or more bands brighter than Bw < 
23.8, R < 22.8, and / < 21.4. Tests with artificial stars indicate < 5% of stars are being 
misclassified as galaxies at the classification limit in all four subfields. 

Regions surrounding saturated stars were removed from the catalog to exclude clus- 
tered spurious objects detected in the wings of the point spread function of bright stars. 
Unclustered contamination reduces the amplitude of the correlation function by (1 — f) 2 
where / is the fractional contamination of the sample. We assume that 1% of bright stars 
are misclassified as galaxies, a rate which is consistent with tests with artificial stars which 
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are a magnitude brighter than the classification limits. While this is a crude assumption, 
misclassified bright stars would only significantly alter the correlation function results if they 
were misclassified at a rate of more than 10%. 

The number of objects classified as stars rapidly decreases at magnitudes fainter than 
R « 22. To estimate the number counts of faint red stars, we assumed stars fainter than 
the classification limits have the same color distribution as stars within a magnitude of the 
classification limits. We then modeled the star counts with a power law given by 

dN 

^— oc 10 am . (1) 
dm 

This relationship was normalized so the counts agreed with the measured star counts at the 
classification limits. For this work, we assumed a = 0.07. While this is only an approxima- 
tion, the increase in the counts is comparable to the observed star counts in the Medium 
Deep Survey (Santiago, Gilmore, & Elson 1996). It is also comparable to the power-law 
index for the counts of 19 < R < 21 red stars in the NDWFS. Tests with a between 0.0 
and 0.15 did not alter the amplitude of the angular correlation function significantly. In 
addition, the amplitude of the correlation function (after correcting for contamination) does 
not significantly change if the star-galaxy classification magnitude limits are decreased. For 
all the correlation function bins, the estimated rate of stellar contamination was less than 
12%. 



3. Photometric redshifts 

Photometric redshifts were determined for all objects with B W) R, and /-band de- 
tections. Model spectra as a function of redshift were generated using PEGASE2 spectral 
synthesis code (Fioc & Rocca-Volmerange 1997). Models with z = Solar metallicity, Miller- 
Scalo initial mass functions, z = ages of 12 Gyr (formation z ~ 4.1), and exponentially 
decreasing star formation rates with e-folding times between 0.6 Gyr and 15 Gyr (r-models) 
for an Q = 0.3 A = 0.7 H$ = 70 km s _1 Mpc^ 1 cosmology 2 were used to estimate galaxy 
colors, /^-corrections and spectral evolution corrections. At z > 3.9, r-models with ages of 
100 Myr were used to model the spectra of high redshift galaxies. The model spectra were 
multiplied by intrinsic dust extinction with E(B — V) = 0.04 and R v = 3.1, comparable to 
estimates for < z < 1 early- type galaxies (Falco et al. 1999; Kauffmann et al. 2003). We 
multiplied the B W) R, and / filter transmission curves with the MOSAIC CCD quantum ef- 
ficiency as a function of wavelength, the mirror reflectivity and a measurement of the KPNO 
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atmospheric extinction to improve the accuracy of the model galaxy colors. Galaxy photom- 
etry was corrected for Galactic dust extinction using the dust maps of Schlegel, Finkbeiner 
& Davis (1998), though it should be noted that the maximum estimate of E(B — V) was only 
0.014 in the four subfields. We used the Vega spectral energy distribution of Hayes (1985) to 
zeropoint the model galaxy colors. Model galaxy colors zeropointed with the Hayes (1985) 
spectrum or the Castelli & Kurucz (1994) model of Vega differ by Byy — R — 0.03 from model 
galaxy colors zeropointed with the frequently used Kurucz (1979) model of Vega. Interpo- 
lation between the r-model was used to fill the color-space occupied by galaxies observed in 
the NDWFS. The uncertainties of the photometric redshifts would be underestimated if the 
color-space was not filled. This is particularly true if color-redshift degeneracies present in 
the data are not reproduced by the models. 

Photometric redshifts were estimated by finding the minimum value of x 2 as a function 
of redshift, spectral type (r), and luminosity. To reduce the CPU time required to evaluate 
the photometric redshifts, the interpolated models were only evaluated when they differed 
from neighboring models sufficiently to significantly alter the photometric redshifts. As the 
model spectra do not account for the observed width of the galaxy locus, we increased the 
photometric uncertainties for the galaxies by 0.05 magnitudes (added in quadrature). This 
results in the photometric redshift code producing la errors which are consistent with the 
observed scatter between the photometric and spectroscopic redshifts in Figure I. 

To further improve the accuracy of the photometric redshifts and their uncertainties, 
the estimated redshift distribution of galaxies as a function of spectral type and apparent 
magnitude was introduced as a prior. This approach, rather than a Bayesian prior (Kodama, 
Bell & Bower 1999; Benitez 2000), was undertaken, as the number of objects as a function of 
r is poorly known at present. The 2dFGRS z = luminosity functions for different spectral 
types (Madgwick et al. 2002), with luminosity evolution given by the r-models, were used to 
produce estimates of the redshift distributions. At M R — 51og/i < —23, where the luminosity 
function of red galaxies is poorly fitted by a Schechter function (Madgwick et al. 2002), 
the luminosity function was approximated by a power-law. To obtain the correspondence 
between the 2dFGRS principal component rj parameter and the r parameter, we fitted 
PEGASE2 models to the 2dFGRS principal component spectra between 3900A and 4100A. 
We did not fit for the entire 2dFGRS spectrum as there may be errors in the 2dFGRS 
continuum calibration (Madgwick et al. 2002). The prior has little effect on the best-fit 
estimates of the photometric redshifts but does significantly alter the photometric redshift 
uncertainties. 

To confirm the reliability of the photometric redshifts, spectral types, and absolute 
magnitudes, simulated data were generated using the r spectral evolution models and the 
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2dFGRS luminosity functions. The simulated data consisted of r < 15 Gyr galaxies with 
the redshift range < z < 5 and luminosity range 0.01 < L* < 100. The simulated object 
photometry was scattered using the estimated uncertainties (including the 0.05 component 
discussed earlier) as a function of apparent magnitude, thus mimicking what would be present 
in the real catalogs. In addition to the simulated data, the accuracy of the photometric 
redshifts was confirmed with spectroscopic redshifts and BwRI photometry for selected 
objects in the NDWFS Bootes, NDWFS Cetus, and Lockman Hole fields. 

A comparison of the photometric and spectroscopic redshifts for —22.5 < Mr — 51og/i < 
— 19.5 red galaxies is shown in right-hand panel of Figure 1. The observed colors of the 
galaxies with spectroscopic redshifts are shown in Figure 2. We discuss the selection criteria 
for the red galaxies in §4. The photometric redshifts have a 5% systematic error and ±5% 
la uncertainty. Only two of the 29 galaxies have photometric redshifts with errors of > 30%. 
One of the outliers has strong [Oil] emission and Balmer absorption. The other is a blended 
object consisting of an emission line galaxy and an M-star. 

Increasing or decreasing the internal dust extinction increases the offset between the 
photometric and spectroscopic redshifts and increases the random scatter. As varying the 
dust extinction does not improve the accuracy of the photometric redshifts, we use the con- 
stant value of E(B — V) — 0.04. Photometric redshifts using the GISSEL01 (Bruzual & 
Chariot 1993; Liu, Chariot & Graham 2000) r-model spectral energy distributions produce 
a slightly larger residual offset and more random scatter, so we use the PEGASE2 models 
throughout the paper. However, there are no large differences between correlation functions 
determined using galaxies with PEGASE2 and GISSEL01 photometric redshifts. Through- 
out the remainder of the paper, photometric redshifts and absolute magnitudes have been 
corrected for the 5% systematic error shown in Figure 1. 



4. The red galaxy sample 

If the PEGASE2 r-models are a good approximation for the spectral evolution of 
z < 0.90 galaxies, the spectral model fits by the photometric redshift code can be used 
to measure the spectral types and luminosity evolution of galaxies. This assumption is con- 
sistent with the B W RI colors of the galaxy locus (Dey et al., in preparation) and comparisons 
of photometric and spectroscopic redshifts discussed in §3. We were therefore able to select 
comparable populations of galaxies at multiple epochs. 

If a galaxy sample is going to be used to measure the evolution of clustering, an evo- 
lutionary sequence of related galaxy populations must be selected at different redshifts. As 
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the clustering of galaxies is known to be a function of absolute magnitude at low redshift 
(e.g., Norberg et al. 2002; Zehavi et al. 2002), accurate photometric redshifts are required 
so objects with similar luminosities can be selected at multiple epochs, enabling unbiased 
studies of the evolution of clustering. Galaxy types with low rates of spectral evolution are 
advantageous as selection criteria relying on rapidly evolving models will be very sensitive 
to errors in the model spectral energy distributions. Red galaxies have accurate photometric 
redshifts and low rates of spectral evolution. At z < 1, the 4000 A break is moving through 
the optical so accurate photometric redshifts can be obtained with a limited number of op- 
tical passbands. The PEGASE2 r ~ 1 Gyr models and observations (J0rgensen et al. 1999; 
Schade et al. 1999; Im et al. 2002) indicate z = 0.9 red galaxies are only « 0.9 magnitudes 
brighter in rest-frame i?-band than z = red galaxies. 

For this work, we selected a sample of red galaxies to be comparable to the early-type 
sample of the 2dFGRS (Norberg et al. 2002). Galaxies fitted with r < 4.5 Gyr were chosen 
so their z = spectra matched the 2dFGRS principal component selection criterion. The 
rest-frame color of the 12 Gyr old r = 4.5 Gyr model with E(B — V) = 0.04 intrinsic dust 
extinction is By/ — R = 1.44. This is approximately the same color as the Sab template 
of Fukugita, Shimasaku, & Ichikawa (1995). The r < 4.5 Gyr selection criterion is ~ 0.15 
magnitudes redder in B w — R than the u* — r* > 1.8 (AB) color cut for SDSS red galaxies 
in Zehavi et al. (2002). A subsample, selected with r < 2.0 Gyr, was chosen to allow 
estimates of clustering as a function of rest-frame color and comparisons with redder samples, 
including z > 0.8 EROs. The rest-frame color of the 12 Gyr old r = 2.0 Gyr model with 
E[B — V) = 0.04 intrinsic dust extinction is B w — R = 1.77, which is only 0.07 magnitudes 
bluer in Bw — R than the 12 Gyr old r = 0.6 Gyr model. The observed colors of the red 
galaxy sample, along with the r-models, are shown in Figure 2. For comparison, the rest- 
frame colors of the Coleman, Wu & Weedman (1980) E and Sbc templates are B^ — R =1.71 
and B\y — R = 1.14. Objects with colors which differ significantly from the models can 
contaminate the sample. To reduce this contamination, only 0.5 < R — I < 1.6 galaxies with 
photometric redshift fits with \ 2 < 3 (including the prior) were included in the final sample. 

The conclusions of this paper rely on the accuracy of the galaxy spectral classifications 
and absolute magnitudes. Significant contamination by late-type spirals or R > 21 faint 
blue galaxies will dramatically decrease angular and spatial two-point correlation functions 
(Efstathiou et al. 1991). To minimize contamination, the photometric redshift range was 
restricted to 0.30 < z < 0.90, where the r-models are not degenerate for B W RI photometry. 
The bulk of the spectra used in Figure 1 have relatively low signal-to-noise ratios making 
it difficult to verify the selection criteria spectroscopically with currently available datasets. 
Instead, we used the simulated data and compared the values of r and Mr used to generate 
the model object with the output r and Mr values from the photometric redshift code. 
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As shown in Figure 3, there is good agreement between the input and output values of r 
and Mr for red galaxies. Few blue galaxies are scattered into the red sample while a small 
fraction of red galaxies are scattered out of the sample. For the remainder of the paper, we 
restrict the magnitude range to R < 23.25, where tests with simulated galaxies indicate r is 
being reliably measured. While some contamination is inevitable, it is reasonable to assume 
that weakly clustered faint blue galaxies (e.g., Efstathiou et al. 1991) are not dominating the 
measured angular correlation function. The final red galaxy sample contains 5325 objects, 
which is 14% of all 39316, R < 23.25 galaxies in the 1.2D sample area. 



5. The correlation function 

We determined the angular correlation function using the Landy & Szalay (1993) esti- 
mator: 

where DD, DR, and RR are the number of galaxy-galaxy, galaxy-random and random- 
random pairs at angular separation 9 ± 56/2. The pair counts were determined in logarith- 
mically spaced bins between 10" and 1°. 

The random objects consist of copies of real galaxies which have had their positions 
changed to mimic objects that are randomly distributed across the sky. This does not result 
in a perfectly uniform surface density of objects across the field due to the completeness 
variations in the 4 subfields. When the "random" objects were distributed across the field, 
the probability of each object being detected in a given subfield was estimated and the 
product of this and the subfield area was used when determining which subfield the object 
would be placed. To decrease the contribution of the random objects to the shot-noise, 100 
random object catalogs were generated and DR and RR were renormalized accordingly. 

The estimator of the correlation function is subject to the integral constraint 

u(e 12 )dQ 1 dQ 2 ~ (3) 



(Groth & Peebles 1977) which results in a systematic underestimate of the clustering. To 
remove this bias, the term 

UJn = ^?J J ^(M^A (4) 

was added to u>(6) where f2 is the survey area. The value of n 2 ujQ, where n is the mean 
number of galaxies per area Q, is the contribution of clustering to the variance of the galaxy 
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number counts (Groth & Peebles 1977; Efstathiou et al. 1991). The angular correlation 
function was assumed to be a power-law given by 

u(0) = u,(l') (J^ 7 (5) 

where 7 is a constant. This is a good approximation of the observed spatial correlation 
function from the 2dFGRS and SDSS surveys on scales of < 10/i -1 Mpc (Norberg et al. 2001; 
Zehavi et al. 2002). For a 7 = 1.87 power-law, the integral constraint for this study was 
approximately 6% of the amplitude of the correlation function at 1'. 

We determined the covariance matrix of 00(8) using the approximation of Eisenstein & 
Zaldarriaga (2001): 

1 f°° 

CM, 0j) = ^pj KP*(K)J (K9 t )J (K9 3 )dK (6) 
where Jo is a Bessel function and P 2 (K) is the angular power spectrum, 

poo 

P 2 (K) = 2tt / w(9)J (K9)9d9. (7) 
Jo 

This approximation is best suited to correlation functions where ui{6) <C 1 and underesti- 
mates the covariance of very strongly clustered objects. For the evaluation of P 2 (K), we 
truncated the power-law form of the correlation function at 2° as the z ~ galaxy corre- 
lation function is ~ on scales of > 20/i~ 1 Mpc (e.g., Maddox, Efstathiou, & Sutherland 
1996; Connolly et al. 2002). However, the power-law fits to the data are only marginally 
affected by the value of P2(K) on large scales. If the angular correlation function bins have 
significant width, Equation 6 is modified to 



CM, 0j) = j2 W j2 W \ \ OO'CM 0')d9d9' (8) 

(D. Eisenstein 2003, private communication) where 9\ and 9 2 are the inner and outer radii 
of the bins. This can be rewritten as the single integral 

4 

Cu{6i > 9j) = (91 -9^(91 -9^ X 
Pl(K) [diMKew) - 9 h MK9 hl )] \9 h2 .h(K9 h2 ) - 9^J X (K9^)\ — . (9) 

The contribution of shot noise to the estimate of the covariance was included by adding the 
reciprocal of the sky surface density of galaxies (per steradian) to P 2 (k). However, the shot 
noise only dominates the covariance on scales of < 1' for the red galaxy sample. 
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The spatial correlation function was obtained using Limber's (1954) equation: 
f°° dN T f°° dN 1 / / f°° dN \ 2 

*[/ ^, M d -^/(l (io, 

where ^ is the redshift distribution, £ is the spatial correlation function and r(9, z, z') is 
the comoving distance between two objects at redshifts z and z' separated by angle 9 on the 
sky. The spatial correlation function was assumed to be a power law given by 

ar,z) = [r/r (z)r (11) 

where 

r (z)=r (0)[l + z]-^-^ (12) 

and e is a constant (Groth & Peebles 1977). Clustering is fixed in physical or comoving 
coordinates ife = 0ore = 7 — 3 respectively. 

We estimated the redshift distribution for Limber's equation by summing the redshift 
likelihood distributions of the individual galaxies in each subsample. As shown in Figure 4, 
the redshift distributions of individual galaxies can not be modeled with Gaussians and 
estimates derived from the photometric redshift code x 2 as a function of redshift must be 
used instead. Model redshift distributions for several subsamples selected by luminosity and 
photometric redshift are shown in Figure 5. 



6. The clustering of red galaxies 

6.1. The angular and spatial correlation functions 

The angular correlation function was determined for red galaxies in a series of photo- 
metric redshift bins between z = 0.30 and z = 0.90, and absolute magnitude bins between 
Mr — 51og/i = —22.5 and Mr — b\ogh = —19.5. All bins are volume limited samples contain- 
ing galaxies brighter than R = 23.25. A power-law of the form u>(6) = u(l')(6 /l') 1-7 was 
fitted to the 9 < 0.25° data with 7 fixed to 1.87, the value for z < 0.15 early-type galaxies 
(Norberg et al. 2002; Zehavi et al. 2002). We use rather than for the power-law 

fits as it depends less on the assumed value of 7. The amplitude of the two-point angular 
correlation functions for these subsamples are summarized in Tables 2 and 3. Angular cor- 
relation functions for —21.5 < M R — 51og/i < 20.5, r < 4.5 Gyr galaxies are also plotted in 
Figure 6. 

A summary of the spatial clustering (parameterized by r ) as a function of spectral type, 
absolute magnitude, and redshift is presented in Tables 2 and 3. The estimates of r in the 
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narrow redshift bins are consistent with the values in the widest redshift bins. While this is 
expected, it is not always the case as the largest redshift bin contains object pairs which are 
not present in the smallest bin. The width of the redshift distribution for the smallest bin 
strongly depends on the uncertainties of the photometric redshifts as these uncertainties are 
comparable to the bin width. In contrast, the shape of the redshift distribution of the largest 
bin depends mostly on the bin width as the uncertainties of the photometric redshifts are 
smaller than the width of the bin. If the uncertainties are systematically underestimated, 
Limber's equation will overestimate the number of close object pairs in the narrowest bin 
and underestimate the value of r . If the uncertainties in the redshift distribution were not 
included, r would vary between 3.9±0.4/i _1 Mpc and 5.5±0.4/i~ 1 Mpc for the narrowest and 
widest redshift bins for —21.0 < M R — 51og/i < —20.0, r < 4.5 Gyr galaxies. The measured 
evolution of r should be independent of the width of the redshift bins and confirming this is 
an extremely useful internal consistency check which should always be applied to correlation 
functions using photometric redshifts. 



6.2. Clustering as a function of absolute magnitude 

Figure 7 and Table 2 present estimates of the spatial clustering of red galaxies as a 
function of evolution corrected absolute magnitude. Absolute magnitude bins containing 
galaxies brighter than M R — 51og/i = —20.5 include the entire 0.30 < z < 0.90 photometric 
redshift range while fainter bins have truncated redshift ranges which are listed in Table 2. 
While redder galaxies are more strongly clustered than bluer galaxies, the striking correlation 
is between r and absolute magnitude. While there is a mild correlation with luminosity at 
M R — 51og/i ~ —20.0, the value of r rapidly increases from r = 6.3/i _1 Mpc at M R — 51og/i = 
—21 to r = 11.2/i -1 Mpc at M R — 51og/i = —22. Similar behavior is seen for at z < 0.15 
in the SSRS2 and 2dFGRS (Willmer, da Costa, & Pellegrini 1998; Norberg et al. 2002), 
and, with lower significance, in CNOC2 at z < 0.4 (Shepherd et al. 2001). Hogg et al. 
(2003) observe similar trends in the SDSS at 0.05 < z < 0.22 by measuring density of galaxy 
neighbors within 8/i^ 1 Mpc spheres as a function of galaxy color and luminosity. Wilson 
(2003) measures r = 4.02 ± 0.22/ i ~ 1 Mpc for 0.2 < z < 1.0, M B < -18.61 red galaxies, but 
her model redshift distribution does not include the uncertainties of her V — I photometric 
redshifts, so her value of r is a lower limit. A summary of previous measurements of red 
galaxy correlation functions is provided in Table 4. 

While the r values of M R - 51og/i > -21.5 galaxies from the NDWFS and 2dFGRS 
show only marginal differences, the r values of luminous galaxies in the NDWFS are higher 
than those of the 2dFGRS with « 2a significance. Clustering evolution would be expected 
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to produce decreasing r values with increasing redshift rather than the opposite trend seen 
in Figure 7. However, as luminous galaxies are more strongly clustered than L* galaxies, 
estimates of their clustering are also more susceptible to cosmic variance. The distribution 
of galaxies on the plane of the sky, which is plotted in Figure 8, clearly shows that luminous 
galaxies are in structures comparable in size to the field-of-view. While it is plausible that 
selection effects could produce the observed structures, stars selected with the same selection 
criteria do not show similar large-scale structure. In addition, in Figure 5, the model redshift 
distribution of the most luminous red galaxies shows evidence of individual structures. We 
therefore assume that the 2a difference between the clustering of luminous galaxies in the 
NDWFS and 2dFGRS is due cosmic variance. Even if cosmic variance were not an issue, it is 
difficult to measure clustering evolution with galaxies in the absolute magnitude range where 
r is strongly correlated with luminosity as small luminosity errors can translate into large 
errors in r . It is therefore preferable to measure clustering evolution with Mr — 51og/i > 
—21.5 galaxies, as we have done in §6.3. 

As the correlation between tq and spectral type is relatively weak in the NDWFS and 
2dFGRS, our measurements of the evolution of clustering are not sensitive to small errors 
in the estimates of the spectral types. The strong correlation between r and absolute 
magnitude for luminous galaxies is a prediction of recent large volume ACDM simulations 
(e.g., Benson et al. 2001). As shown in Figure 7, a 141 3 /i~ 3 Mpc 3 ACDM simulation with 
galaxy selection criteria similar to the r < 4.5 Gyr sample (Benson et al. 2001, A. Benson 
2002, private communication) is a good approximation of the clustering of M R ~ — 21 red 
galaxies. 



6.3. Evolution of the spatial correlation function 

The evolution of spatial clustering galaxies was studied with —21.5 < M R — 51og/i < 
—20.5 red galaxies. Luminous galaxies were excluded as their spatial clustering is strongly 
correlated with absolute magnitude and small redshift or spectral evolution errors could 
produce large changes in the measured spatial clustering. The faint limit was chosen to 
allow the same range of absolute magnitudes to be studied over a broad redshift range. 

As shown in Figure 9 and Table 3, no significant evolution of r (comoving) occurs 
over the redshift range studied. Two models of the clustering of red galaxies, derived from 
the GIF ACDM simulations (Jenkins et al. 1998), are plotted in Figure 9 and provide good 
approximations to the measured clustering from the NDWFS, SDSS and 2dFGRS. The 
Kauffmann et al. (1999) model predicts the clustering of early-type galaxies with stellar 
masses of > 3 x 1O 1O M while the Benson (2002, private communication) simulation models 
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the clustering of B - R > 1.24, -20.5 < M R - 51og/i < -21.5 galaxies. If the evolution of 
the underlying dark matter distribution at z < 1 is well described by the linear or quasi- 
linear growth of density perturbations, the bias of red galaxies must be rapidly evolving with 
redshift. 

The evolution of r was empirically measured by estimating the clustering evolution 
parameter e (from Equation 12). The 0.30 < z < 0.50 and 0.70 < z < 0.90 bins were 
assumed to be independent estimates of r at the median redshifts of their model redshift 
distributions. The 2dFGRS and SDSS clustering estimates were also included to provide 
additional constraints on e. The selection criteria for the NDWFS r < 4.5 Gyr sample 
allow direct comparison with the 2dFGRS early-type sample of Norberg et al. (2002). The 
present-day colors of the NDWFS r < 4.5 Gyr sample are only ~ 0.15 magnitudes redder in 
B w - R than the SDSS u* - r* > 1.8 sample of Zehavi et al. (2002). The best-fit estimates 
of e are summarized in Table 5 and plotted in Figure 10. Models with e = are rejected with 
> 2.5cr confidence when constraints from the 2dFGRS are included with the NDWFS data. 
The evolution of tq is consistent with clustering fixed in comoving coordinates (e = 7 — 3) 
at z < 0.90. 



7. Discussion 

The clustering of red galaxies undergoes little or no evolution from z ~ 0.9 to the 
present-day. At z > 0.9, the only published clustering study using comparable template 
selection criteria is that of Firth et al. (2002) which used galaxies fitted with the evolving E 
and Sbc templates of the HYPERZ photometric redshift code (Bolzonella, Miralles & Pello 
2000). Their estimate of r = 7.0 ± 1.9/i _1 Mpc, with 7 assumed to be 1.8, is consistent with 
tq remaining fixed in comoving coordinates to z ~ 1.5. However, to differentiate between 
currently plausible ACDM models larger survey areas, on the order of the complete NDWFS 
(?s 18D°), will be required. 

Further constraints on the clustering of 1 < z < 2 galaxies are available from samples 
of QSOs (Croom et al. 2001) and R — K > 5 extremely red objects (EROs; Daddi et al. 
2001; Firth et al. 2002; Roche et al. 2002). The QSO spatial correlation function, measured 
in redshift space (s replacing r ), marginally increases from s = 5.28 ±o;|g /i -1 Mpc with 
7 = 1.72±°;| at z < 0.95 to s = 6.93 ±{f 4 h^Mpc with 7 = 1.64±°;|> at 2.10 < z < 2.90 
(Croom et al. 2001). However, as pointed out by Croom et al. (2001), interpretation of QSO 
clustering relies on poorly constrained models of QSOs lifetimes and host populations as a 
function of redshift. 
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EROs are easier to relate to low redshift populations than QSOs as the majority appear 
to be the progenitors of early-type galaxies with the remainder being dusty starbursts (Dey 
et al. 1999; Liu et al. 2000; Moriondo et al. 2000). EROs could therefore extend the redshift 
range of early-type spatial correlation function evolution estimates to z ~ 1.5. With the 
exception of H < 20.5 EROs in the Las Campanas Infrared Survey (Firth et al. 2002), all 
EROs clustering studies measure r > 10/i _1 Mpc (Daddi et al. 2001; Firth et al. 2002; Roche 
et al. 2002). This is comparable to the clustering of the most luminous red and early-type 
galaxies at z < 1. 

If K > 19 EROs correspond to < L* galaxies (Daddi et al. 2001; Roche et al. 2002), 
the ERO clustering results are difficult to reconcile with z < 1 studies and ACDM theory. 
The clustering of red L* galaxies in the NDWFS can be increased to r ~ 9/i _1 Mpc only by 
increasing the photometric redshift la uncertainties to ~ 40% of the photometric redshifts 
or increasing the stellar contamination to ~ 40%. These scenarios are inconsistent with the 
redshift comparisons in Figure 1 and the lack of a stellar locus in Figure 2. The discrepancy 
between the ERO values and z < 1 studies could be due to errors in the ERO model 
redshift distributions used to deproject the angular correlation function. Firth et al. (2002), 
who measure r = 7.7 ± 2.4/i _1 Mpc for H < 20.5 EROs, have few z > 1.5 objects while 
the Daddi et al. (2001) and Roche et al. (2002) models have a significant fraction of EROs 
at z pa 2. Improved constraints on ERO clustering should be provided with model redshift 
distributions constrained with photometric redshifts which have been verified with spectro- 
scopic samples. The clustering of EROs, as measured from the NDWFS, will be described 
in a future paper (Brown et al., in preparation) using a B W RIK dataset. 



8. Summary 

We have used the NOAO Deep Wide-Field Survey to measure the clustering of 0.30 < 
z < 0.90 red galaxies. The wide-field and B W RI bands allow large galaxy samples to be 
selected as a function of spectral type and absolute magnitude using photometric redshifts. 
PEGASE2 spectral evolution models with exponentially decreasing star formation rates have 
been used to select an evolutionary sequence of related galaxies as a function of redshift. The 
red sample, with present-day rest-frame colors of Bw — R> 1.44, was chosen to allow direct 
comparison with the low redshift early-type sample from the 2dFGRS. The clustering of red 
galaxies is strongly correlated with luminosity, with r increasing from 4.4 ± 0.4/i _1 Mpc at 
M R - 51og/i « -20.0 to 11.2 ± 1.0/i _1 Mpc at M R - 51og/i « -22.0. Clustering evolution 
measurements with samples where the distribution of spectral types and luminosities are a 
function of redshift will be dominated by selection effects. The strength of r (comoving) 
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as a function of absolute magnitude in our sample is comparable to estimates at z < 0.15 
from the 2dFGRS, with differences at high luminosity being attributable to structures of sizes 
comparable to the field-of-view. No significant evolution of r was detected in comparisons of 
the NDWFS with the 2dFGRS and SDSS. For 0.30 < z < 0.90, -21.5 < M R -h\ogh < -20.5, 
Bw — R > 1.44 galaxies, the largest sample studied, the value of tq is 6.3 ± 0.5/i _1 Mpc with 
7 fixed at 1.87. The strong clustering and lack of detectable evolution appears consistent 
with recent ACDM models where the bias undergoes rapid evolution and r undergoes little 
evolution at z < 2. 
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Spectroscopic Redshift Spectroscopic Redshift 

Fig. 1. — A comparison of the photometric and spectroscopic redshifts of —22.5 < Mr — 
51og/i < —19.5 red galaxies (r < 4.5 Gyr) with photometric redshifts in the range 0.30 < 
z < 0.90. Dashed horizontal lines mark the upper and lower limits of the photometric 
redshift range. Dotted diagonal lines show the median offset and ±lcr uncertainty of the 
photometric redshifts. Simulated galaxies, generated using the r-models with photometric 
noise added, are shown in the lefthand panel. On the right are real galaxies with spectroscopic 
redshifts from the NED database and unpublished spectra. The real galaxies show a 5% 
systematic offset between the photometric and spectroscopic redshifts. With the inclusion of 
0.05 magnitudes of intrinsic scatter, the real data have la uncertainty estimates which are 
consistent with the measured scatter between the photometric and spectroscopic redshifts. 
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Fig. 2. — The Vega colors of the r < 4.5 Gyr, 0.30 < z < 0.90 NDWFS Bootes galaxy sample 
(left panel) and the r-models (right panel). Bold points in the left hand panel show the colors 
of galaxies with spectroscopic redshifts (including Cetus and Lockman Hole galaxies) which 
were used to test the accuracy of the photometric redshifts. Bold points on the model tracks 
mark the predicted colors of z = 0, 1, 2 and 3 galaxies for each of the models. For clarity, 
only five of the r-model tracks are plotted. The r-models occupy the same region of color 
space as the NDWFS red galaxies. Single burst models (not shown) have different colors as 
a function of redshift. A more detailed description of the colors of NDWFS galaxies and the 
accuracy of galaxy stellar evolution models will be provided by Dey et al. (in preparation). 
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Input t (Gyr) Input M R -5log/i 

Fig. 3. — A comparison of the input and measured spectral types and absolute magnitudes for 
simulated galaxies with photometric redshifts in the range 0.30 < z < 0.90. The r < 4.5 Gyr 
selection criterion is shown with the dotted line in the left-hand panel. There is good 
agreement between the input and output r and Mr values for red galaxies. Few blue galaxies 
are scattered into the red galaxy sample. Some structure can be observed in the r output 
values due to the interpolated models only being evaluated when their colors and redshifts 
differ significantly from the other models. 
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Redshift 

Fig. 4. — The redshift likelihood distribution produced by the photometric redshift code. 
NDWFS J142518.9 + 353545, which has a spectroscopic redshift of 0.705, has a photometric 
redshift of 0.63 (shown by the arrow) and a spectral type of r = 0.64 Gyr. Unlike Brunner 
et al. (2000) and Brown, Boyle & Webster (2001), we can not adequately model the redshift 
likelihood distribution with a Gaussian. 
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Fig. 5. — The redshift distributions used to evaluate the spatial correlation function of the 
r < 4.5 Gyr samples as a function of luminosity (left panel) and redshift (right panel). The 
redshift distributions of adjoining photometric redshift bins have significant overlaps and 
the resulting estimates of r$ for these bins cannot be considered entirely independent. The 
redshift distribution of the most luminous red galaxies shows evidence of individual structures 
in redshift space. The redshift distributions of the r < 2.0 Gyr subsamples (which are not 
shown) have comparable shapes but different normalizations. 
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Fig. 6. — The angular correlation functions of r < 4.5 Gyr, —21.5 < Mr — hlogh < —20.5 
galaxies selected from the NDWFS. The first three panels show the angular correlation func- 
tions for the narrowest photometric redshift bins while the last panel shows the angular 
correlation function for the 0.30 < z < 0.90 bin. On small angular scales shot noise domi- 
nates the uncertainties while on larger scales the contribution of clustering to the covariance 
dominates. Power-law fits to the data are shown along with ±lcr errors (dotted lines). For 
reference, 0.1° corresponds to a transverse comoving distance of 3/t _1 Mpc at z ~ 0.7. 
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Fig. 7. — The spatial clustering of red galaxies as a function of absolute magnitude. For 
clarity, only the uncertainties for the r < 4.5 Gyr data are shown and the r < 2.0 Gyr data 
has been shifted slightly to the right. While datapoints brighter than M« — 51og/i = —20.5 
include galaxies over the entire 0.30 < z < 0.90 photometric redshift range, fainter bins have 
a truncated redshift range. The z < 0.15 2dFGRS data points, which have uncertainties of 
~ 0.5/i~ 1 Mpc, are over-plotted using the assumption that Bj — R = 1.25 for red galaxies. 
The clustering of red galaxies is a strong function of luminosity and a weaker function of 
spectral type. A ACDM simulation of the clustering B — R > 1.24 galaxies (Benson et al. 
2001, A. Benson 2002, private communication) reproduces the clustering of L* galaxies but 
the 141 3 /i~ 3 Mpc 3 simulation volume is insufficient to measure the clustering of luminous 
galaxies. 
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Fig. 8. — The distribution of red galaxies on the plane of the sky for several abso- 
lute magnitude bins. The bins are (from left to right) —20.5 < Mr — 51og/i < —19.5, 
-21.5 <M R - 51og/i < -20.5 and -22.5 <M R - 51og/i < -21.5. The greyscale of the dots 
is related to the photometric redshift, with dark dots at low redshift and light dots at high 
redshift. Rectangular voids in the galaxy distribution are the boundaries between the ND- 
WFS subfields and regions surrounding bright stars. The distribution of luminous galaxies 
clearly shows evidence for structures comparable to the area of the subfields, suggesting that 
surveys of this size are susceptible to cosmic variance. 
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Fig. 9. — The spatial correlation function parameter r$ as a function of redshift for the ND- 
WFS, 2dFGRS, and SDSS red galaxy samples. For clarity, the r < 2.0 Gyr data points have 
been shifted slightly to the right. The horizontal bars indicate the photometric (NDWFS) 
and spectroscopic (2dFGRS, SDSS) redshift ranges. The clustering of r < 4.5 Gyr galaxies 
(which were selected to allow direct comparison with the 2dFGRS) is well approximated by 
ACDM simulations of the clustering of early-type and B — R > 1.24 galaxies (Kauffmann 
et al. 1999; Benson et al. 2001, A. Benson 2002, private communication). From the plot, it 
is clear the value of ro (in comoving coordinates) undergoes little or no evolution over the 
redshift range observed. 
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Fig. 10. — The evolution of r (comoving). The NDWFS data are shown with the com- 
parable low redshift data from the 2dFGRS (left panel) and SDSS (right panel). Models 
with clustering fixed in physical coordinates (e = 0) are rejected with 2.5a confidence by the 
NDWFS and 2dFGRS while models with little or no evolution (e ~ 7 — 3) provide acceptable 
fits to the observed clustering. 
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Table 1. The four subflelds of the first NDWFS public data release. 



Subfield name 


R.A. 


Decl. 


Delivered FWHM (") a 


Integration timc(hours) 


50% completeness limit b 




(J2000.0) 


(J2000.0) 


B w 


R 


I 


B w 


R 




% 


R 


I 


NDWFS J1426+3531 


14 26 00.8 


+35 31 32 


1.2 


1.6 


1.3 


2.1 


1.7 


3.0 


26.8 


24.8 


24.7 


NDWFS J1426+3456 


14 26 01.4 


+34 56 32 


1.3 


1.2 


1.0 


2.3 


1.2 


2.1 


26.5 


25.3 


24.9 


NDWFS J1428+3531 


14 28 52.8 


+35 31 39 


1.5 


1.3 


0.7 


2.3 


1.7 


2.2 


26.4 


25.0 


25.6 


NDWFS J1428+3456 


14 28 52.2 


+34 56 39 


1.6 


1.3 


1.3 


2.3 


1.7 


2.6 


26.2 


25.0 


23.6 



a Measured using the SExtractor FWHM values of bright unsaturated stars. 

b Determined with SExtractor using artificial stellar objects inserted into copies of the real data. 



Table 2. The correlation function of red galaxies in the NDWFS as a function of absolute 

magnitude. 



Selection 


Photo, z range 


Absolute magnitude 


Apparent magnitude 


Galaxies 


Contamination 


w(l') 


Median z 


roO^Mpc) 


T < 4.5 Gyr 


0.30-0.90 


-22.50 < M R < -21.50 


18.18 < R < 22.43 


660 


0.009 


0.41 ±0.07 


0.64 


11.2 ± 1.0 


t < 4.5 Gyr 


0.30-0.90 


-22.00 < M R < -21.00 


18.30 < R < 23.01 


1677 


0.005 


0.24 ±0.03 


0.61 


8.5 ±0.6 


t < 4.5 Gyr 


0.30-0.90 


-21.50 < M R < -20.50 


18.76 < R < 23.25 


2651 


0.023 


0.14 ±0.02 


0.60 


6.3 ±0.5 


t < 4.5 Gyr 


0.30-0.75 


-21.00 < M R < -20.00 


19.17 < R < 23.17 


2429 


0.009 


0.12 ±0.02 


0.55 


5.0 ±0.4 


t < 4.5 Gyr 


0.30-0.65 


-20.50 < M R < -19.50 


19.72 < R < 23.08 


1756 


0.025 


0.12 ±0.02 


0.51 


4.4 ±0.4 


t < 2.0 Gyr 


0.30-0.90 


-22.50 < M R < -21.50 


18.18 < R < 22.43 


488 


0.011 


0.52 ±0.10 


0.66 


12.4 ± 1.2 


t < 2.0 Gyr 


0.30-0.90 


-22.00 < M R < -21.00 


18.30 < R < 23.01 


1122 


0.006 


0.29 ±0.05 


0.64 


9.2 ±0.8 


r < 2.0 Gyr 


0.30-0.90 


-21.50 < M R < -20.50 


18.76 < R < 23.25 


1592 


0.037 


0.20 ±0.03 


0.62 


7.5 ±0.6 


t < 2.0 Gyr 


0.30-0.75 


-21.00 < M R < -20.00 


19.25 < R < 23.17 


1337 


0.016 


0.16 ±0.04 


0.57 


5.7 ±0.7 


t < 2.0 Gyr 


0.30-0.65 


-20.50 < M R < -19.50 


19.87 < R < 23.08 


929 


0.043 


0.14 ±0.04 


0.53 


4.7 ±0.7 



Table 3. The correlation function of red galaxies in the NDWFS as a function of 

photometric redshift. 



Selection Photo, z range Absolute magnitude Apparent magnitude Galaxies Contamination ^(1') Median z ro(h Mpc) 



t < 4.5 


Gyr 


0.30-0.50 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


21.25 


853 


0.003 





.25 ±0.05 


0.42 


5.3 ±0.6 


t < 4.5 


Gyr 


0.50-0.70 


-21 


.50 


< 


Mr 


< 


-20.50 


20.05 


< 


R < 


22.50 


1023 


0.002 





.24 ±0.05 


0.60 


5.5 ±0.6 


t < 4.5 


Gyr 


0.70-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


20.99 


< 


R < 


23.25 


775 


0.071 





.20 ±0.05 


0.85 


6.7 ±0.8 


t < 4.5 


Gyr 


0.30-0.70 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


22.50 


1876 


0.003 





.19 ±0.03 


0.52 


5.9 ±0.5 


t < 4.5 


Gyr 


0.50-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


20.05 


< 


R < 


23.25 


1798 


0.032 





.17 ±0.03 


0.69 


6.3 ±0.5 


t < 4.5 


Gyr 


0.30-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


23.25 


2651 


0.023 





.14 ±0.02 


0.60 


6.3 ±0.5 


t < 2.0 


Gyr 


0.30-0.50 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


21.25 


435 


0.005 





.31 ±0.08 


0.44 


5.5 ±0.8 


t < 2.0 


Gyr 


0.50-0.70 


-21 


.50 


< 


Mr 


< 


-20.50 


20.12 


< 


R < 


22.50 


665 


0.004 





.29 ±0.07 


0.61 


5.9 ±0.8 


t < 2.0 


Gyr 


0.70-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


21.20 


< 


R < 


23.25 


492 


0.111 





.31 ±0.07 


0.85 


8.4 ± 1.0 


t < 2.0 


Gyr 


0.30-0.70 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


22.50 


1100 


0.004 





.20 ±0.04 


0.55 


5.8 ±0.7 


t < 2.0 


Gyr 


0.50-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


20.12 


< 


R < 


23.25 


1157 


0.050 





.22 ±0.04 


0.69 


7.1 ±0.7 


t < 2.0 


Gyr 


0.30-0.90 


-21 


.50 


< 


Mr 


< 


-20.50 


18.76 


< 


R < 


23.25 


1592 


0.037 





.20 ±0.03 


0.62 


7.5 ±0.6 



Table 4. A summary of several studies of early-type and red galaxy correlation functions. 



Survey a ' 6 


Redshift Range 


Galaxies 


Magnitude Range 


Selection 


rrj(/i J Mpc comoving) c ' d 


7 C 


NDWFS 


0.30 < z < 0.90 


2651 


-21.5 < Mr < -20.5 


t < 4.5 Gyr 


6.3 ± 0.5 


1.87 


NDWFS 


0.30 < z < 0.90 


1592 


-21.5 < Mr < -20.5 


t < 2.0 Gyr 


7.5 ± 0.6 


1.87 


Perseus-Pisces 


2 < 0.04 


278 


M Zw < -19.5 


Morphology 


8.35 ±0.75 


2 n5+ 010 


SSRS2 


z < 0.020 


395 


M B < -19.4 


Morphology 


5.10 ±0.38 


1.91 ±0.26 


SSRS2 


2 < 0.027 


418 


M B < -20.0 


Morphology 


5.27 ±0.46 


1.86 ±0.29 


SSRS2 


2 < 0.033 


372 


M B < -20.5 


Morphology 


5.73 ±0.56 


2.30 ± 0.46 


SSRS2 


2 < 0.040 


272 


M B < -20.9 


Morphology 


8.60 ± 1.44 


2.45 ±0.71 


APM 


2 £ 0.1 


336 


-22 < M B < -15 


Morphology 


5.9 ±0.7 


1.85 ±0.13 


2dFGRS 


0.016 < 2 < 0.071 


1909 


-18.5 < Mb r < -17.5 


Red SED (rj < -1.4) 


8.33 ± 1.82 


1.87 ±0.23 


2dFGRS 


0.020 < 2 < 0.086 


3717 


-19.0 < Mb r < -18.0 


Red SED (ri < -1.4) 


6.28 ± 1.46 


1.98 ±0.11 


2dFGRS 


0.031 < 2 < 0.125 


10,135 


-20.0 < Mb , < -19.0 


Red SED (r? < -1.4) 


5.71 ±0.57 


1.87 ±0.09 


2dFGRS 


0.048 < 2 < 0.150 


6434 


-21.0 < Mb j < -20.0 


Red SED (r? < -1.4) 


6.10 ±0.72 


1.80 ±0.12 


2dFGRS 


0.072 < 2 < 0.150 


686 


-22.0 < Mb j < -21.0 


Red SED (r) < -1.4) 


9.74 ± 1.16 


1.95 ±0.37 


SDSS 


0.019 < 2 < 0.13 


19,603 


14.5 < r*'< 17.6 


Rest-frame u* — r* > 1.8 


6.78 ±0.23 


1.86 ±0.03 



PDF 


2 < 0.5 


22,359 


Bj < 22.5 


Redder than CWW Sbc 


6.0 ±0.3 


1.9 


K20 


0.5 < 2 < 2.0 


400 


K < 19.2 


R- K > 5 


12 ±3 


1.8 


CNOC2 


0.120 < 2 < 0.270 


248 


R c < 21.5, M R < -20 


Red CWW templates 


5.35 ±0.20 


2.05 ±0.08 


CNOC2 


0.270 < 2 < 0.382 


234 


R c < 21.5, M R < -20 


Red CWW templates 


6.55 ± 1.16 


2.05 ±0.08 


CNOC2 


0.382 < 2 < 0.510 


238 


R c < 21.5, Mr < -20 


Red CWW templates 


6.99 ±0.57 


2.10 ±0.11 


CNOC2 


0.12 < 2 < 0.40 


254 


-20.00 < Mr < -19.25 


Red CWW templates 


5.82 ±0.81 


1.85 ±0.08 


CNOC2 


0.12 < 2 < 0.40 


276 


-20.65 < Mr < -20.00 


Red CWW templates 


5.40 ±0.19 


1.89 ±0.06 


CNOC2 


0.12 < 2 < 0.40 


278 


-22.52 < Mr < -20.65 


Red CWW templates 


6.71 ±0.53 


2.14 ±0.02 


UH8K 


0.20 < 2 < 0.90 


3382 


M B < -18.61 


Redder than CWW Sbc 


4.02 ±0.22 


1.8 


LCIRS 


0.3 <2 <0.8 


272 


H < 20.5 


Evolving E & Sbc 


7.0 ± 1.6 


1.8 


LCIRS 


0.8 < 2 < 1.5 


355 


H < 20.5 


Evolving E & Sbc 


7.0 ± 1.9 


1.8 


LCIRS 


0.7 < 2 < 1.5 


337 


H < 20.0 


R-H >4 


11.1 ± 2.0 


1.8 


LCIRS 


0.7 < 2 < 1.5 


312 


H < 20.5 


R- H > 4 


7.7 ±2.4 


1.8 


ELAIS N2 


0.5 < 2 < 2.0 


166 


K < 21.0 


R- K > 5 


10-13 


1.8 



a Perseus-Pisces (Guzzo ct al. 1997), SSRS2 (Willmer, da Costa, & Pellegrini 1998), APM (Loveday et al. 1995), 2dFGRS (Norberg et al. 2002), 
SDSS (Zehavi et al. 2002), PDF (Brown, Boyle & Webster 2001), K20 (Daddi et al. 2001), CNOC2 (Shepherd et al. 2001), UH8K (Wilson 2003), 
LCIRS (Firth et al. 2002), and ELAIS N2 (Roche et al. 2002). 

b For clarity the surveys have been split into 3 groups; NDWFS, z < 0.2, and z > 0.2. 

c For other studies, uncertainties are as published and may not include the effect of the covariance on the uncertainty estimates. 
d The parameter rJJ/1.73 is used instead of r for CNOC2 (Shepherd ct al. 2001). 
°Where the value of 7 was fixed, the value is given without an error estimate. 



Table 5. Constraints on e (as defined in Equation 12) provided by the NDWFS, SDSS, 

and 2dFGRS. 



NDWFS Sample z < 0.2 Sample a r (0)(/i _1 Mpc comoving) e 



T 


< 2.0 Gyr 




none 




3.1±i;3b 




T 


< 4.5 Gyr 




none 




3.8±\i h 


-2.9±}f 


T 


< 4.5 Gyr 


SDSS u* - r* 


> 1.8, -22.0 < M r * < 


-19.0 


6.9±8;| 


-o-5±8:S 


T 


< 4.5 Gyr 


2dFGRS r, < 


-1.4, -20.0 < M B j < 


-19.0 


5 S+ - 7 
°- : ' =l= 0.6 


-i-6±8:l 



a The median redshifts of the SDSS and 2dFGRS samples are assumed to be 0.1. We use the uncertainties 
of Norberg et al. (2002) and Zchavi ct al. (2002) as published. 

b As the NDWFS red galaxy sample photometric redshift range is 0.30 < z < 0.90, the value of ro(0) from 
the NDWFS alone is an extrapolation, which is strongly dependent on the best-fit value of e. 



